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PREDICTION  OF  EXPLOSIVELY  DRIVEN 
RELATIVE  DISPLACEMENTS  IN  ROCKS 

Scott  Blouin 


INTRODUCTION 

Deep-based  defense  systems,  such  as  a  reserve  ballis¬ 
tic  missile  force  or  command  control  and  communica¬ 
tions  facilities,  buried  in  rock  at  depths  ot  halt  a  kilo¬ 
meter  01  more,  hold  promise  as  relatively  invulnciable 
future  links  in  United  States  strategic  defenses.  One 
of  the  kev  uncertainties  in  deep  basing,  however,  is 
the  possible  susceptibility  ot  these  systems  to  relative 
displacements  of  rock,  or  block  motions,  which  occur 
well  beyond  the  crater  caused  by  surface  or  shallow- 
buried  nuclear  explosions.  Considerable  eisperimental 
evidence  of  rock  displacement  has  been  gathered 
following  both  high  explosive  and  underground  nudeai 
detonations. 

Experimental  studies  of  the  feasibility  of  suing 
missile  silos  in  rock  during  the  Hard  Rock  Silo  Program 
-conducted  by  the  Air  Force  and  the  Defense  Atomic 
Support  Agency  in  the  late  60'sand  early  70's-demon- 
strated  that  a  principal  threat  to  the  silos  was  from 
explosion-induced  displacement  of  large  rock  masses 
or  blocks  along  pre-existing  planes  of  weakness  in  the 
in  situ  rock.  These  relative  motions  were  observed 
along  bedding  planes,  joints,  joint  sets,  and  interfaces 
between  different  rock  types,  in  both  igneous  and 
sedimentary  rocks.  In  instances  where  they  intersected 
concrete  and  steel-lined  structures,  the  structures 
were  sheared  in  direct  proportion  to  the  displacement 
in  the  free-field. 

The  high  explosive  test  senes  which  drove  these 
displacements  is  described  by  Blouin  (1980).  The 
lest  technique  (having  the  acronym  DIHEST,  for 
Direct  Induced  High  Explosive  Simulation  Technique) 
employed  a  buried  vertical  array  of  high  explosives, 
simultaneously  de’r-r'ited,  to  generate  ground  motions 


thought  to  be  similar  to  those  created  by  the  crater¬ 
ing  action  of  a  nuclear  surface  burst.  The  DIHEST 
shots  all  produced  elliptically  shaped  craters  with 
their  long  axes  in  line  with  the  explosive  arrays.  Rel¬ 
ative  displacements  were  observed  to  nearly  three 
crater  radii  from  the  arrays.  They  were  also  observed 
at  depths  below  the  bottom  of  the  craters,  though 
no  attempt  was  made  to  determine  the  maximum 
extent  of  these  displacements  at  depth.  It  is  impor¬ 
tant  to  note  that  of  the  multitude  of  relative  displace- 
m -nts  documented  during  this  program,  perhaps  only 
one  would  have  been  observed  during  a  casual  post- 
shot  examination.  The  rest  were  discovered  following 
the  extensive  debris  clearance  and  excavation  needed 
to  recover  survey  points  on  the  test  bed  surface  or  to 
document  structural  damage.  The  magnitudes  of  the 
relative  displacements  observed  during  this  test  series 
ranged  from  less  than  1  cm  to  4  m. 

The  recent  underground  nuclear  shot  MIGHTY 
EPIC  conclusively  demonstrated  that  explosively 
driven  relative  displacements  of  considerable  magni¬ 
tude  occur  at  great  depths.  Short  and  Kennedy  (1978) 
and  Townsend  ( 1 976,  1 977  a,b,cj  document  at  least 
seven  instances  of  relative  displacement  along  faults, 
bedding  planes,  and  an  interface  between  adjoining 
rock  types.  The  magnitudes  of  these  displacements 
ranged  trom  0.4  m  to  more  than  1 .7  m.  Several  test 
structures  were  severely  damaged  by  these  displace¬ 
ments. 

Clearly,  the  viability  of  any  deep-based  defense 
system  must  partially  rest  on  a  credible  demonstra¬ 
tion  that  the  system  will  not  be  vulnerable  iO  these 
types  of  relative  displacement.  Design  must  follow 
one  or  a  combination  of  three  approaches.  First,  the 
system  mav  be  located  below  the  maximum  depth 
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.  to  which  block  motion  will  be  propagated  by  the  ex¬ 
pected  threat  (surface  or  shallow-buried  bursts). 
Second,  the  system  may  be  designed  to  "absorb”  the 
anticipated  relative  displacements.  Third,  the  rock 
discontinuities  might  be  reinforced  using  rock  bolts, 
grout  or  some  other  means  to  reduce  or  eliminate  the 
anticipated  displacements.  All  three  approaches, 
however,  are  dependent  on  predictions  of  the  extent, 
magnitude,  and  direction  of  the  relative  displacements, 
or  on  the  shear  stresses  which  generate  them.  In  addi¬ 
tion  design  of  critical  communications  links,  access 
tunnels  and  shafts,  etc.  will  also  depend  on  a  credible 
block-motion  prediction  technique.  In  short,  such 
predictions  will  be  critical  to  feasibility  studies  and 
subsequently  to  nearly  all  aspects  of  the  design  of  any 
deep-based  system. 

This  report  summarizes  an  initial  attempt  to  for¬ 
mulate  a  semi-empirical  procedure  for  predicting  the 
maximum  extent,  direction  and  magnitude  of  rela¬ 
tive  displacements  resulting  from  an  underground 
nuclear  detonation  in  the  Nevada  Test  Site  tuff. 
Specifically,  such  predictions  are  derived  for  every 
.fault  documentation  station  of  the  underground  nu¬ 
clear  shot  DIABLO  HAWK,  detonated  in  September 
1978.  This  effort  is  envisioned  as  a  principal  logical 
step  toward  the  ultimate  goal  of  providing  credible 
predictions  of  the  relative  displacements  at  depth 
which  result  from  surface  or  shallow-buried  nuclear 
explosions. 

The  prediction  procedures  described  herein  are 
based  on  relative  displacement  data  from  the  DIHEST 
HE  (high  explosive)  events  and  the  MIGHTY  EPIC 
underground  nuclear  shot.  All  displacements  used  in 
the  analyses  arc  of  the  "driven”  variety;  that  is,  the 
direction  of  slip  is  in  the  general  direction  of  the  max¬ 
imum  dynamic  shear  stress  generated  by  the  explosion. 
So  called  “triggered”  displacements,  where  the  direc¬ 
tion  of  motion  is  in  the  direction  of  the  maximum 
in  situ  shear  stress,  are  not  treated  here  (for  a  further 
discussion  of  triggered  motions,  see  Bache  and  Lambert 
1976).  Analyses  of  both  the  relevant  DIHEST  data 
and  MIGHTY  EPIC  data  are  presented,  along  with  the 
many  assumptions  and  approximations  involved  in 
the  formulation  of  the  prediction  procedures.  Finally, 
details  of  the  DIABLO  HAWK  prediction  procedures 
are  summarized  and  the  individual  predictions  at  each 
fault  documentation  station  are  presented. 

The  predictions  are  made  in  two  parts.  First, 
probable  slip  locations  and  directions  are  identified 
through  an  analysis  procedure  which  superimposes 
an  approximation  of  the  dynamic,  explosion-generated 
stresses  on  the  in  situ  stresses.  Once  probable  slip 
locations  are  identified,  the  magnitudes  of  the  dis¬ 
placements  expected  at  these  locations  are  determined 
through  a  semi-empirical  conservation  of  energy  re¬ 
lationship. 


DIHEST  ANALYSIS 

The  DIHEST  series  produced  significant  block 
motions  in  five  out  of  seven  events.  Though  no 
attempt  was  made  to  monitor  displacements  beneath 
the  explosive  arrays,  the  near-surface  data  are  valuable 
in  that  a  fairly  complete  definition  of  the  extent  and 
geometry  of  the  block  movements  could  be  made 
from  the  combined  test  bed  and  structural  surveys, 
a  definition  which  is  much  more  complete  than  is 
possible  in  the  underground  environment.  In  addi¬ 
tion,  the  large  number  of  ground-motion  instrument 
locations  provided  good  ground-motion  data  on  both 
sides  of  the  slip  planes,  whereas  the  underground  data 
are  more  sparse,  tend  to  be  of  poorer  quality  and  are 
inconclusive  (see,  for  example,  Short  and  Kennedy 
1978).  An  example  of  the  DIHEST  ground  motion 
data  which  are  suggestive  of  the  mechanisms  control¬ 
ling  the  magnitude  of  those  displacements  is  shown 
in  Figure  1.  A  comparison  of  horizontal  velocity 
time  histories  is  shown  for  near-surface  gage  locations 
on  the  STARMET  event.  The  gage  at  the  9.1 -m  range 
from  the  explosive  array  was  located  within  a  large 
block  of  rock  which  was  upthrust  along  a  joint  or 
joint  set  dipping  toward  the  array  at  an  angle  of  67°. 
The  gage  at  the  18.3-m  range  was  located  beyond  this 
slip  surface.  It  is  obvious  from  the  time  histories  that 
relative  motion  began  at,  or  soon  after,  passage  of  the 
initial  velocity  peak-which  should  correspond  closely 
to  passage  of  the  peak  dynamic  stress.  By  the  time 
the  18.3-m  gage  had  come  to  rest  (at  about  0.1  s) 
the  9.1 -m  gage  was  still  moving  outward  at  nearly 
peak  velocity.  This  velocity  was  then  gradually 
slowed  at  a  nearly  constant  rate,  suggesting  applica¬ 
tion  of  a  constant  restraining  force  to  the  sliding 
block.  Total  permanent  horizontal  displacement  at 
the  18.3-m  gage  was  negligible,  while  that  at  the  close- 
in  gage  was  0.51  m.  Both  displacements  agree  closely 
with  post-shot  survey  data. 

In  this  case  the  data  comparison  suggests  that 
once  relative  displacement  begins,  the  thrust  block 
may  be  considered  a  free  body  sliding  up  the  slip  sur¬ 
face  at  an  initial  velocity  equal  to  the  vector  sum  of 
the  peak  horizontal  and  vertical  velocities.  The  con¬ 
stant  restraining  force  is  supplied  by  gravity  acting 
on  the  block,  plus  a  friction  force  between  the  block 
and  slip  plane  resisting  the  upward  slide.  Such  a 
block  model  is  shown  schematically  in  Figure  2. 

The  kinetic  energy  of  the  block  is  thus  expended  in 
the  form  of  work  against  the  forces  of  gravity  and 
friction.  The  work  done  against  gravity  can  be  ex¬ 
pressed  as  a  change  in  potential  energy  of  the  block, 
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Figure  2.  Schema*  ot  1','HCSr  block 
model. 


given  by 

5sin^  (|) 

where  m  is  the  truss  of  the  block,  g  the  acceleration 
of  gravity,  and  h  sin£  the  vertical  component  ot  total 
displacement  along  the  slip  plane,  whetc  (?  is  the  dip 
angle  of  the  clip  surface.  The  work  dune  against  fric¬ 
tion  equals  the  friction  force  multiplied  by  the  total 
displacement  along  the  slip  plane,  which  is  g;ven  by 

ii,  -  y  mg  h  cos  3  (21 

where  y  is  the  coefficient  of  sliding  friction  between 
the  block  and  slip  plane.  The  normal  Knee  on  the 
slip  plane,  mg  cos  0,  is  assumed  to  be  due  only  to  the 
weight  of  the  bluck.  Summing  eq  1  and  2  and  setting 
them  equal  tu  the  initial  kinetic  energy  of  the  block 
yields 

Vmtv^  =  mg  f>  sinp  +  y  mg  u  cos,3  (3) 

where  vn  is  the  vector  sum  of  the  peak  horizontal  and 
vertical  velocities.  Rearranging  gives  an  expression  for 


the  total  displacement  fi  of 

v} 

6  =  u 

2 y(y  costf  +  sin  (!) 

Table  1  summarizes  the  most  dramatic  relative 
Ois.pl.  ements  observed  in  the  DIHEST  senes.  Motion 
generally  occurred  along  joints  or  joint  sets  striking 
approximately  paral'el  to  the  explosive  array  and 
dipping  toward  it.  PLANE  WAVE  II  was  an  excep¬ 
tion  in  that  motion  occurred  along  bedding  planes  and 
lithologic  interfaces  having  essentially  no  dip.  HANDEC 
II  was  also  an  exception,  as  motion  occurred  along  a 
joint  dipping  away  from  the  explosive  array.  The 
last  column  of  Table  i  shows  the  relative  displace¬ 
ments  calculated  using  eq  4.  The  resultant  peak 
velocities  used  in  the  calculations  and  listed  in  the 
table  were  calculated  using  the  average  peak  horizon 
tai  and  vertical  velocities  at  the  range  where  each 
slip  plane  intersected  the  testbed  surface.  The  co¬ 
efficient  of  friction  was  assumed  at  0.5,  a  typical 
value  for  many  faults. 

f  or  all  events  but  DATEX  II,  eq  4  overpredicts 
the  total  displacement.  This  would  be  expected  for 
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Table  1.  DIHEST  relative  displacement  summary. 


Eve  nt 

Rock  type 

Seismic  velocity  (tn/s) 

DIHEST 

yield 

(kg) 

(tons) 

Max  ret 
dlspl  (m) 
Range  (m) 

Dip  angle  0 

A  vg.  peak  particle  velocity 
(mis) 

Horizontal  Vertical  Resultant 

Displ  from 
eg  4 

R  =  .S 
(m) 

PLANEWAVE  II 

<0.1 

Intcrbedded  sedimentary 

1905 

17.4 

2.0* 

- 

>2.0 

0.4** 

1220  -  3050 

(2.1) 

S0° 

STARMET 

1.7 

Granite 

1995 

12.8 

4.6 

4.6 

6.5 

1.9 

3650  -  5 ISO 

(2.2) 

67° 

DATEX  II 

4.0 

Tonalite 

37,190 

33.5 

4.6* 

2.3* 

5.1* 

1.6 

2130  -  3660 

(41) 

—  20° 

HANDEC  II 

0.8 

Tonalite 

41,725 

36.6 

6.1 

2.9 

6.8 

2.5** 

2130  -  3660 

(46) 

-  30° 

ROCKTEST  II 

0.3 

Tonalite 

106,120 

68.6 

2.9 

1.5 

3.3 

0.7 

2130  -  3660 

(117) 

20° 

*  Estimated  value. 

**  Equation  applied,  though  recognized  that  assumptions  governing  its  derivation  do  not  hold  in  this  use. 


PLANEWAVE  II  and  HANDEC  II,  where  the  respec¬ 
tive  horizontal  and  negative  dip  of  the  slip  plane  geo¬ 
metry  tends  to  restrict  slip  more  than  the  geometry 
of  the  model  used  to  derive  the  equation.  The  over¬ 
prediction  on  STARMET  is  minimal;  the  agreement 
between  calculated  and  actual  displacement  is  ex¬ 
cellent.  A  factor  contributing  to  the  overprediction 
of  displacement  for  both  ROCKTEST  II  and  HANDEC 
II  was  that  during  displacement  a  considerable  air 
overpressure  had  been  applied  to  the  testbeds  to 
simulate  the  concurrent  nuclear  airblast.  In  the  case 
of  HANDEC  II,  it  is  esimated  (Blouin  1980)  that 
air  overpressure  on  the  testbed  was  close  to  70  bars 
at  the  arrival  of  the  compression  wave  from  the 
DIHEST  explosion  and  that  the  overpressure  had 
dropped  to  approximately  IS  bars  by  the  completion 
of  relative  displacement.  The  ROCKTEST  II  over¬ 
pressure  is  classified,  but  it  too  was  substantial.  The 
combined  effects  of  the  airblast  and  airblast-induced 
ground  motions  on  the  relative  displacements  are 
complicated  and  beyond  the  scope  of  this  work. 
Qualitatively  at  least,  it  appears  likely  that  they 
tended  to  reduce  the  magnitude  of  slip,  though  not 
to  the  extent  the  first  order  approximations  would 
indicate. 

A  closer  look  at  the  STARMET  ground  motion 
data  reveals  that  the  good  agreement  between  cal¬ 
culated  and  measured  relative  displacements  Is  some¬ 
what  misleading,  but  at  the  same  time  it  Is  indicative 


of  a  mechanism  which  would  explain  the  large  under- 
prediction  of  displacement  on  DATEX  ll-where, 
unfortunately,  no  ground  motion  measurements  were 
taken  in  the  vicinity  of  the  displacements.  Figure  3 
shows  the  vertical  velocity  time  history  from  hole  15 
within  the  STARMET  thrust  block,  the  companion 
to  the  horizontal  time  history  of  Figure  1 .  Starting 
at  approximately  0.1  s  (just  beyond  the  velocity 
peak)  and  lasting  for  more  than  0.5  s,  a  gradual 
steady  deceleration  of  the  block  occurred.  As  indi¬ 
cated,  the  slope  of  this  portion  of  the  time  history  is 
-Ifi,  implying  that  gravity  alone  was  responsible  for 
tf,c  deceleration  of  the  block.  The  displacement 
hodograph  of  Figure  4  was  constructed  from  the 
horizontal  and  vertical  displacement  time  histories. 

It  shows  that  the  thrust  block  moved  upward  and 
outward  along  the  slip  plane,  which  dipped  67®  to¬ 
ward  the  explosive  array,  reaching  a  peak  displace¬ 
ment  at  this  point  of  2.0  m.  The  block  then  slid  back 
down  the  boundary  joint,  coming  to  rest  at  a  final 
displacement  of  1 .4  m.  Using  the  peak  2.5-  and  5.3- 
m/s  horizontal  and  vertical  velocities  from  the  time 
histories  in  Figures  1  and  3  in  eq  4  gives  a  value  of 
displacement  of  1 .9  m  for  a  friction  coefficient  of 
zero  and  a  value  of  1 .6  m  for  a  coefficient  of  0.5. 

The  first  value  is  In  excellent  agreement  with  the 
2.0-m  actual  peak  displacement  while  the  second 
value  agrees  well  with  the  1 ,4-m  permanent  displace¬ 
ment,  However,  it  Is  obvious  that  the  latter  agreement 
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Figure  3.  STARMET  thrust  block-hole  IS  vertical  velocity. 


Figure  4.  Section  schematic  -hole  15  displacement 
trajectory. 

is  simply  fortuitous  because  the  block  slid  back  down 
the  slip  plane  and  happened  to  finally  come  to  rest 
at  a  1.4-m  displacement.  Thus  it  appears  that  the 
block  was  actually  in  spall,  driven  upward  and  out¬ 
ward  by  the  explosion  along  the  direction  of  the  slip 
surface,  but  with  little  or  no  normal  force  on  the 
joint  during  the  outward  transition. 

If  this  assumption  is  applied  to  the  DATEX  II 
maximum  displacement,  eq  4  then  indicates  a  3.9-m 
peak  displacement  when  using  zero  for  the  friction 
coefficient.  This  is  very  close  to  the  actual  perma¬ 
nent  displacement  of  4.0  m.  Furthermore,  for  a  co¬ 
efficient  of  friction  of  0.5,  the  minimum  slip  surface 
dip  which  would  allow  a  block  to  slide  back  down 
the  slip  surface  under  the  force  of  gravity  is  26.6°. 
Thus,  because  of  its  low  dip  angle,  the  DATEX  II 
block  would  not  have  slid  back  down  the  boundary 
joint  as  the  STARMET  block  did,  making  its  perma¬ 
nent  displacement  and  peak  displacement  essentially 
equal. 

In  conclusion,  it  appears  that  the  conservation  of 
block  energy,  using  the  peak  resultant  velocity  to 
calculate  the  initial  kinetic  energy  of  the  block,  is  a 


valid  approach  for  predicting  the  displacement  of 
blocks  which  intersect  the  free  surface.  Due  to  the 
explosion  and/or  ground  motion  geometry,  the 
frictional  resistance  to  sliding  may  be  negligible,  and 
in  some  instances  the  blocks  appear  to  be  in  a  state 
of  spall.  Use  of  eq  4  with  a  coefficient  of  friction  of 
zero  serves  as  an  upper  bound  to  all  the  DIHEST  rel¬ 
ative  displacements.  For  STARMET  and  DATEX  II, 
the  upper  bound  is  a  good  fit  to  the  actual  peak  dis¬ 
placements.  For  HANDEC  II  and  ROCKTEST  II, 
where  the  DIHEST  test  was  combined  with  a  nuclear 
airblast  simulation,  the  interaction  of  the  air  over¬ 
pressure  and  of  the  resulting  airblast-induced  ground 
motions  with  the  DIHEST-induced  relative  displace¬ 
ments  appear  to  have  somewhat  restrained  the  block 
motion. 


MIGHTY  EPIC  ANALYSIS 

Geological  setting  and  relative  displacement 
documentation 

The  MIGHTY  EPIC  underground  nuclear  detona¬ 
tion  produced  six  documented  instances  of  relative 
displacement  along  faults  and  bedding  planes  in  the 
Rainier  Mesa  tuff.  In  most  instances  both  the  dip  and 
strike  components  of  slip  were  determined.  In  addi¬ 
tion,  relative  movement  was  detected  along  the  inter¬ 
face  between  the  tuff  and  an  underlying  quartzite, 
but  since  this  motion  was  within  or  adjacent  to  a 
rubble  zone  which  is  not  close  to  the  location  of  the 
DIABLO  HAWK  event,  It  was  not  incorporated  into 
the  MIGHTY  EPIC  relative  displacement  analysis. 

Only  a  brief  description  of  the  MIGHTY  EPIC/ 
DIABLO  HAWK  geologic  setting  and  relative  displace¬ 
ments  Is  included  here.  For  a  full  description  see  the 
summary  by  Short  and  Kennedy  (1978).  Both 
MIGHTY  EPIC  and  DIABLO  HAWK  were  sited  at 
the  Nevada  Test  Site  in  the  U  12n.  10  tunnel  system 
beneath  Rainier  Mesa.  The  tunnel  system,  shown  in 
plan  view  In  Figure  5,  is  approximately  400  m  be¬ 
neath  the  Mesa  surface  and  extends  horizontally  Into 
the  shot  area  from  the  portal  In  the  side  of  the  mesa. 
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Figure  5.  MIGHTY  I'/’IC  runnel  system  and  muior  laults  itrcm  Shutl  and  Kama 
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figure  6.  Schematic  section  tut  t  .quart?  ite  interface  (from  Short  and 
Kennedy  1 97K). 
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The  two  shot  locations  (or  working  points)  were  about 
1 52  m  apart,  lying  in  the  Main,  or  Line  of  Site  (LOS) 
drift.  MIGHTY  EPIC  was  fired  in  1976  and  DIABLO 
HAWK  in  1 978.  The  A,  B  and  C  drifts  contained  ex¬ 
perimental  structures  of  various  configurations  and 
hardnesses.  Additional  tunneling  and  experiments 
"were  added  prior  to  DIABLO  HAWK.  The  interface 
drift  provided  access  to  the  tuff/quartzite  interface 
which  dipped  south  beneath  the  drift  and  the  working 
point.  The  top  of  the  competent  quartzite,  as  shown 
schematically  in  Figure  6,  lay  approximately  19  m 
below  the  floor  of  the  drift  at  hole  MH3, 46  m  below 
the  floor  at  hole  MH1 ,  and  79  m  below  the  working 
point. 

The  material  properties  of  the  tuff  in  the  area  of 
the  drifts  containing  the  structures  are  reported  by 
Butters  (1976).  The  tuff  had  an  average  porosity  of 
31%,  a  density  of  1.95  g/cm3,  and  an  unconfincd 
compressive  strength  of  195  bars.  Ultrasonic  labora¬ 
tory  p  and  s  wave  velocities  averaged  3140  and  1 580 
m/s  respectively.  Tabulated  data  presented  by  Short 
and  Kennedy  (1978)  indicate  that  average  tuff  prop¬ 
erties  in  other  areas  generally  vary  less  than  10%  from 
these  values. 

Figure  5  shows  a  plan  view  of  the  major  faults  in 
the  area  of  the  MIGHTY  EPIC  works.  Fault  numbers 


3  through  8  are  throughgoing  and  can  be  extrapolated 
to  faults  or  lineations  on  the  surface  of  the  mesa.  All 
these  major  faults  are  normal  faults,  and  most  strike 
north  to  northwest,  dip  steeply  toward  the  west,  and 
exhibit  small  displacements.  The  faults  are  very  tight, 
with  the  fault  zones  generally  less  than  0.01  m  thick. 
Short  and  Kennedy  (1978)  report  that  the  faults 
themselves,  as  well  as  the  surrounding  rock,  have  very 
low  pcrmeability-usually  well  below  a  millidarcy.  In 
addition  to  the  faults,  there  are  also  numerous  bedding 
planes  within  the  tuff.  These  have  a  shallow  dip  (on 
the  order  of  1 5°)  toward  the  south-southeast. 

In  situ  stresses  beneath  Rainier  Mesa  appear  to  vary 
considerably  in  magnitude  and  orientation.  Figure  7 
is  a  plot  of  maximum  and  minimum  horizontal  stress 
from  seven  locations  beneath  the  mesa,  obtained  using 
an  overcoring  technique  (Ellis  1976).  According  to 
Short  and  Kennedy  (1978),  the  maximum  horizontal 
stress  averaged  70  bars,  the  minimum  horizontal  stress 
averaged  33  bars,  and  the  vertical  stress  averaged  63 
bars,  which  agrees  with  the  overburden  stress  at  the 
depth  of  the  measurements.  The  inclination  of  the 
actual  principal  stresses  varies  from  site  to  site.  The 
maximum  horizontal  stresses  are  oriented  toward  the 
northeast  and  the  minimum  horizontal  stresses  toward 
the  northwest. 


Figure  7.  Secondary  principal  stresses  In  horizontal  plane  for 
seven  locations  In  Rainier  Mesa  (from  Ellis  1976). 
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For  all  MIGHTY  EPIC  and  DIABLO  HAWK  stress 
calculations  in  this  study,  the  approximation  of  in 
situ  stresses  shown  in  Figure  8  was  used.  The  principal 
stresses  were  assumed  to  have  vertical  and  horizontal 
orientations.  A  vertical  stress  of  80  bars  was  assumed. 
This  was  approximately  equal  to  the  overburden  stress 
at  the  depth  of  the  shots.  The  maximum  horizontal 
principal  stress  (80  bars)  was  assumed  equal  to  the 
vertical  stress  and  oriented  N  45°  E.  The  minimum 
horizontal  principal  stress  (40  bars)  was  assumed 
equal  to  half  the  maximum  and  was  oriented  N  45°  W. 

The  relative  displacements  documented  for  MIGHTY 
EPIC  are  shown  in  the  plan  view  of  Figure  9  and  are 
listed  in  Table  2  (both  after  Kipp  and  Kennedy  1978). 
All  displacements;  except  the  fault  near  the  end  of 
the  Interface  drift,  were  similar  to  those  observed 
during  the  DIHEST  events  in  that  relative  displacement 
was  always  relative  to  the  foot  wall,  as  shown  sche¬ 
matically  in  Figure  10.  It  appears  likely  that  upward 
vertical  motion  along  the  tuff/quartz  interface  beneath 
the  interface  fault  resulted  in  displacement  in  the 
normal  mode  along  this  fault. 

•  Relative  displacementsalong  the  interface,  as  deter¬ 
mined  by  a  post-test  magnetometer  survey,  were  ap¬ 
proximately  2.1  m  outward  from  the  working  point 
and  upward  at  the  61 -m  range  of  hole  MH1  (see  Fig. 

6),  dropping  to  about  0.8  m  at  the  1 22-m  range  of 
hole  MH3.  Both  the  remainder  of  the  fault  displace¬ 
ments  and  the  bedding  plane  displacement  between 
tunnel  bed  subunits  3BC  and  3D  had  reverse  dip  slip. 
The  magnitude  of  displacement  of  the  small  fault 
between  faults  5  and  6  which  intersected  the  SRI  mini¬ 
structure  could  not  be  determined.  Structural  defor¬ 
mation  indicated  that  the  fault  moved  with  reverse 
dip  slip  and  left  lateral  strike  components.  The  strike 


Figure  8.  Shot  level  approximation 
of  principal  In  situ  stresses. 


slip  components  of  the  3BC-3D  bedding  plane  slip 
and  the  fault  near  fault  6  in  the  LOS  drift  could  not 
be  determined. 

The  displacement  indicated  in  Figure  9  at  the  inter¬ 
section  of  fault  6  with  the  Bypass  drift  is  somewhat 
confusing.  According  to  Kipp  and  Kennedy  (1978), 
the  fault  at  this  location  is  actually  a  zone  of  multiple 
discontinuities  rather  than  a  well  defined  slip  surface. 
Relative  displacement  was  measured  only  in  the  right 
rib  (or  north  side)  of  the  drift.  This  apparently  con¬ 
sisted  of  an  0.8-m  compression  across  the  fault  zone 
and  a  0.6-m  increase  in  drift  diameter.  There  is  a 
possibility  that  displacement  in  this  area  was  rotational. 
Because  of  these  irregularities,  this  motion  was  not 
listed  in  Table  2,  nor  is  ii  included  in  the  MIGHTY 
EPIC  analysis. 

Shear  stress  analysis 

Analysis  of  the  MIGHTY  EPIC  relative  displacements 
was  separated  into  two  parts.  First,  the  approximate 
peak  shear  stresses  acting  at  the  location  of  each  dis¬ 
placement  were  analyzed  in  an  attempt  to  define  the 
properties  governing  the  in  situ  shear  strengths  of  the 
faults.  Second,  the  magnitudes  of  displacement  were 
analyzed  to  devise  a  semi-empirical  model  for  predic¬ 
ting  fault  displacements  at  depth.  The  first  part  of 
the  analysis  resulted  in  a  method  for  predicting  both 
the  location  and  direction  of  relative  displacements, 
the  second  resulted  in  a  method  for  predicting  their 
magnitude. 

The  prediction  of  relative  displacement  locations 
and  directions  is  based  on  the  calculation  of  the  peak 
shear  and  normal  stresses  acting  on  a  particular  fault, 
at  a  particular  location  on  that  fault.  The  peak  shear 
stress  is  then  compared  to  the  critical  shear  stress 
given  by  the  Mohr-Coulomb  failure  criterion  as 

Tc  =  T0  +  °n  Un^  (5) 

where  rc  is  the  critical  shear  stress,  which  is  equiva¬ 
lent  to  the  shear  strength  of  the  fault  under  the  con¬ 
ditions  of  in  situ  stress  and  dynamic  loading,  rQ  is  a 
measure  of  the  cohesion  and/or  mechanical  inter¬ 
locking  of  the  fault  under  zero  normal  stress,  on  is 

Table  2.  Observed  MIGHTY  EPIC  relative  displacements. 

Ditploctmeni  (ml 

Slip  loco  lion  Strike  flip  Dip  flip  Total 


r au It  5  B  drift 

0.40 

0.24 

0.46 

f  ault  5  Bypaas  drift 

0.4$ 

0.49 

0.67 

1  ault  near  fault  6-LOS  drift 

- 

1.66 

>1.6S 

Ml  _/3D  bedding  plane 

'  - 

0.91 

>0.91 

fault  through  SRI  mlnittructurt 

>0 

>0 

>0 

fault  through  Interface  drift 

•  0.76 

ai0.46 

ftOJt 

8 


■ 


legend 

•  SBC  «CK1  und*r  3D  3' 

•  •  Down  J.6* 

M  No  movtm«nt 


r iyurc  9.  Block  motion  resulting  from  the  MltjH  CY  EPIC,  event  (from  Short  and 
Kennedy  1 9  78). 


Normal  Dirtciion  Rovorso  Direction 


f  igure  W.  Schematic  section  view-normal  and  reverse  faults. 


and 


the  peak  normal  stress  acting  on  the  fault,  and  0  is 
the  friction  angle  of  th'  fault.  If  the  peak  shear  stress 
exceeds  the  critical  shear  stress  rc  it  is  assumed  that 
motion  along  the  fault  plane  will  result.  The  direction 
of  slip  would  be  expected  to  follow  the  direction  of 
the  peak  shear  stress  vector.  The  magnitude  of  peak 
shear  is  given  by 


where  ts  and  rd  are  the  peak  shear  stresses  in  the 
direction  of  the  fault  strike  and  dip,  respectively,  and 
t,  is  the  resultant  of  these  two  shear  vectors-with  a 
magnitude  equal  to  their  Pythagorean  sum  and  a 
direction  governed  by  their  signs. 

The  expressions  for  the  peak  normal  stress  and  peak 
dip  and  strike  shear  stresses  were  derived  by  using  the 
schematic  of  the  slip  plane  shown  in  Figure  1 1  and 
notation  from  Timoshenko  and  Goodier  (1970).  The 
slip  plane,  represented  by  the  plane  abc,  is  oriented 
with  respect  to  the  orthogonal  axes  x,  y  and  /,  along 
which  the  peak  principal  stresses  act  as  shown.  For 
the  geometry  of  the  MIGHTY  EPIC  and  DIABLO 
HAWK  experiments,  with  the  fault  documentation 
locations  at  the  same  level  as  the  shot  point  and  the 
principal  in  situ  stresses  oriented  horizontally  and 
vertically  as  depicted  in  Figure  8,  the  peak  combined 
principal  stresses  will  also  act  horizontally  and  verti¬ 
cally.  Thus,  the  angle  6  is  taken  as  the  angle  between 
the  slip  plane  strike  and  the  maximum  peak  horizontal 
principal  stress  op  which  acts  in  the  x  direction.  The 
minimum  peak  horizontal  principal  stress  acts  in 
the  y  direction  and  the  vertical  principal  stress  a v 
acts  in  the  z  direction.  The  angle  0  is  the  slip  plane 
dip.  Line-on  is  perpendicular  to  the  slip  plane  at 
point  n.  Let  X,  Y  and  Z  be  components  of  stress  on 
plane  abc  acting  in  the  principal  stress  orientation  x, 
y  and  z,  respectively,  and  let  the  area  of  abc  be  equal 
to  A. 

To  satisfy  equilibrium,  the  sum  of  forces  in  each 
of  the  x,  y  and  z  directions  must  equal  zero.  The  sum 
of  forces  in  the  x  direction  is  expressed  as 

=  XA  -  A  cos  [aon)  op  =  0  (7) 

where  A  cos  (aon)  equals  the  area  of  plane  obc.  Solv¬ 
ing  for  X  yields 

X  =  op  cos  (aon)  (8) 

and  a  similar  procedure  for  forces  in  they  and  z 
directions  yields 

Y=  om  cos  (bon)  (9) 


Z  =  ov  cos  (con)  (10) 

Using  the  stress  components  X,  Y  and  Z,  the  shear 
stress  in  the  direction  of  the  slip  plane  strike  r,  the 
shear  stress  in  the  direction  of  the  slip  plane  dip  rd, 
and  the  normal  stress  on  the  slip  plane  op  can  be 
obtained  as  shown  in  the  horizontal  and  vertical 
sections  through  the  slip  plane  in  Figure  12.  The 
strike  shear,  defined  as  positive  when  acting  on  the 
slip  plane  in  the  negative  x  direction,  is  given  by 

r  =  X  cos  0  -  Y  sinC  (11) 

Combining  with  cq  8  and  9  yields 

T  -  a  cos  (aon)  cos 0  -  o_  cos  (bon)  sin  0 

\  r>  '  fTl 

(12) 

The  dip  shear,  defined  as  positive  when  acting  upward 
on  the  slip  plane,  is  given  by 

r.  ~  X  <ut0  cos  P*  Y  cos  -  Z  sin  0  (13) 

a 

Combining  with  eq  8-10  yields 

Trf  =  o  cos  (aon)  sin  0  cos 0  + 
o  cos  (bon)  cos  0  cos0  -  a  cos  (con)  sin  0  (14) 

The  normal  stress  on  the  slip  plane,  defined  as  posi¬ 
tive  when  acting  toward  the  origin,  is  given  by 

o  -  X  cos  (aon)  +  Y  cos  (bon)  +  Z  cos  (con) 
"  (15) 

which  combined  with  eq  8-10  yields 

o=o  cos 2  (aon)  +  om  cos2  (bon)  + 
n  p  ’  m 

cos2  (con)  (16) 

The  cosines  of  aon,  bon  and  con  must  be  expressed 
in  terms  of  0  and  0  From  Figure  1 1 : 

cos  (aon )  =  ™  (17) 

where 

on  =  o/7?sin0  (18) 

and 

om  =  oo  sin  6  (19) 
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figure  1 1.  Geometry  of  slip  plane  abc. 


Substituting  cq  18  and  19  into  cq  1  7  yields 

cos  ( aon )  =  sin  0  sin  p  (20) 

Similarly,  it  can  be  shown  that 

cos  (bon)  =  cos  0  sin  p  (21 ) 

and  by  inspection, 

cos  (con)  =  cos p  (22) 

Substituting  eq  20-22  into  cq  12,  14  and  1 6  and  manip 
Dialing  ihem  yields  the  following  expressions  for  the 
shear  and  normal  stresses  on  the  slip  plane  as  functions 
of  the  peak  piincipal  stresses  and  the  fault  dip  and 
strike: 

r  =  0P'°m  (sin  20  sin  (5)  (23) 


Slip  plane -horizontal  section 


z 


Slip  plane-vertical  section 

Figure  1 2.  Sections  through  slip  plane 
abc. 


lop  s,n 


cos 


2  0  -  o  ) -ii 
2 

(24) 


on  =  op  sin2  0  sin2  p  +  om  cos  2  U  sin2  p  + 

ov  cos2  p  (25) 

In  order  to  apply  the  above  equabons  at  a  particu¬ 
lar  location  along  any  fault  or  bedding  plane,  the 
magnitude  and  orientation  of  the  peak  principal 
stresses  must  be  established.  These  are  obtained  by 
superimposing  an  approximation  of  the  critical  princi¬ 
pal  dynamic  stresses  at  the  location  of  interest  on  the 
in  situ  stress  approximation  shown  in  Figure  8.  The 
critical  dynamic  slate  of  stress  is  assumed  to  occur 
shortly  after  passage  of  the  stress  front.  The  peak 
principal  dynamic  stress  is  radial  in  compression  and 
is  approximated  by 


0,-pc  vt  (26) 

where  p  is  the  mass  density  of  the  rock,  c  the  com¬ 
pressive  wave  speed,  and  v,  the  peak  radial  particle 
velocity  at  the  location  of  interest  on  the  slip  plane. 
Figure  13  is  a  logarithmic  plot  of  peax  particle  veloc¬ 
ity  and  stress,  scaled  by  the  cube  root  of  yield,  used 
to  determine  peak  radial  stresses  for  both  MIGHTY 
EPIC  and  DIABLO  HAWK  (from  Short  and  Kennedy 
1978).  This  plot  is  a  fit  to  past  data  in  this  rock  type 
and  it  proved  to  be  a  good  fit  to  the  MIGHTY  EPIC 
data.  A  compression  wave  speed  of  2590  m/s  and 
density  of  1 .96  g/cnr'  were  used  to  convert  the 
peak  particle  velocity  to  stress.  The  minor  dynamic 
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Figure  13.  MIGHTY  FPIC-DIABLO  HAWK 
peak  radial  stress  and  particle  velocity  approx¬ 
imation  (after  Short  and  Kennedy  1978). 


pi  inti  pal  sIicsjcs  die  in  the  tangential  direction.  Un 
fortunately ,  according  to  Bass  (pers.  comm,),  there 
arc  no  credible  tangential  stress  data  upon  which  to 
base  approximations  of  the  critical  minor  principal 
stresses.  Thus,  the  elastic  solution  for  a  spherically 
expanding  stress  wave  from  a  spherical  cavity  was  used 
as  a  guide  in  picking  the  critical  minor  stresses. 

Selbcrg  (1951)  and  Rinehart  (1975)  give  solutions 
for  a  Heaviside  pressure  puise  (nondecaying  step  in¬ 
put),  applied  to  a  sphere  of  radius  z0  in  an  infinite 
elastic  medium  having  a  compression  wave  speed  c 
and  Poisson's  ratio  of  0.25.  Graphs  of  radial  stress, 
tangential  stress  and  maximum  shear  stress  for  ranges 
of  fg,  2 and  infinity  are  shown  in  Figure  19  (Irom 
Selbcrg  1951).  F  or  comparison,  the  curves  arc  nor¬ 
malised  by  multiplying  the  amplitudes  b\  r;r  and 
by  defining  time  t  as 


where  t  is  real  time. 

At  the  wave  front,  both  radial  and  tangential  stresses 
are  compressive.  This  is  followed  by  a  rapid  icversal 
of  tangential  stress  into  tension.  At  a  range  of  2 rQ, 
which  most  nearly  corresponds  to  the  ranges  of  inter¬ 
est  m  the  underground  tests,  the  maximum  shear 
stress  remains  nearly  constant  during  the  course  of 


Tim#  Parameter  (f'| 

b  Tantjtmioi  Stress 


Tii*#  Paramal#'  (»') 

C-  Mflumuffi  Srt«ar  Strtss 


Figure  14.  \ormahrcd  stresses 
from  a  spherically  expanding 
wove  in  an  elastic  medium  (after 
.S elberg  1932  and  Rinehart  1973). 


this  reversal ,  while  at  smaller  ranges  it  increases  and 
at  greaicr  ranges  .1  decreases.  Based  on  these  obser¬ 
vations,  a  logical  approximation  of  the  critical  dy namic 
stress  slate  (i.e.,  that  state  producing  the  highest  com¬ 
ponents  of  shear  stress!  is  radial  stress  at  its  peak 
value.  This  peak  value  of  radial  stress  occurs  shortly 
behind  the  initial  arrival  ut  the  compression  wave  and, 
with  tangential  stress  at  zero,  would  represent  the 
rapid  reversal  in  tangential  stress.  T ensile  stresses  are 
not  allowed  since  it  is  believed  that  large  masses  of  in 
situ  rock  can  support  little  or  no  tensile  stress  because 
of  the  many  zones  of  weakness  which  are  always  pre¬ 
sent.  In  actuality  this  assumption  of  critical  stresses 
is  probably  conservative  since  the  radial  stress  should 
begin  decreasing  by  the  time  tangential  stress  drops 
to  zero. 

With  the  critical  dynamic  stress  state  now  defined, 
it  must  be  superimposed  on  the  in  situ  stresses  to  ob¬ 
tain  the  magnitude  and  orientation  of  the  principal 
stresses  acting  on  the  slip  plane  at  the  point  t  f  inter¬ 
est.  For  both  MIGHTY  EPIC  and  DIABLO  HAWK, 
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the  displacement  documentation  locations  arc  at  the 
same  elevation  as  the  working  point.  Therefore,  the 
dsnamic  radial  stress  has  no  component  in  the  vertical 
direction  and  will  interact  only  with  the  horizontal 
components  of  the  principal  stress  approximation 
shown  in  Figure  8.  A  schematic  plan  view  of  this  in¬ 
teraction  is  shown  in  Figure  IS.  Angle  q  is  defined  as 
the  angle  between  the  direction  of  propagation  of  the 
radial  stress  of  from  the  working  point  and  the  normal 
to  the  plane  on  which  the  principal  stress  acts.  It  is 
measured  counterclockwise.  Mohr  circles  are  used  to 
depict  the  superposition  of  stresses  In  Figure  16.  The 
first  portion  of  the  figure  shows  the  dynamic  critical 
principal  stresses  of  and  the  assumed  critical  tangen¬ 
tial  stress  of  zero.  Diameter  ob  leprcsents  the  princi¬ 
pal  planes  on  which  the  in  situ  stresses  act  The  nor¬ 
mal  component  ot  d>  nanne  stress  acting  on  the  major 
in  situ  principal  plane  a  is  gisen  by 


Figure  1 5.  Interaction  between  dynamic  and  in  situ 
principal  stresses. 


Figure  16.  Mohr's  circle  representat’on-summation  of  horizontal  stresses. 


°rp  ~2 —  (I  +  cos  2a)  (28) 

and  the  normal  component  of  dynamic  stress  acting 
on  the  minor  in  situ  principal  plane  arm  is  given  by 

°rrn  ~  ~ —  0  -COs2<*)  (^ 

The  shear  component  of  dynamic  stress  acting  on  the 
major  in  situ  principal  plane  7rp  is  given  by 

TrP  =~2~  sin2a  (3°) 

and  the  dynamic  shear  component  acting  on  the  minor 
plane  rrm  is  given  by 

r  =  - _ —  sin2a  (31) 

rm  2 


is  obtained  by  subtracting  this  expression.  Figure  1 7 
is  a  schematic  view  of  the  relationship  among  the  com¬ 
bined  principal  stresses,  the  original  in  situ  stresses 
and  the  dynamic  stress. 

The  combined  principal  stresses  can  then  be  used 
in  eq  23-25  to  define  the  peak  normal  and  shear 
stresses  on  any  potential  slip  plane.  The  angle  0  be¬ 
tween  the  major  principal  stress  and  the  strike  of  the 
slip  plane  is  determined  by  finding  the  orientation  of 
the  principal  stresses  using  eq  32.  The  peak  strike  and 
dip  shear  stress  are  combined  using  eq  6  to  obtain  the 
resultant  peak  shear  stress  on  the  slip  plane.  Peak 
resultant  shear  stresses  were  thus  obtained  at  each 
MIGHTY  EPIC  relative  displacement  location  listed 
in  Table  2.  Since  failure  occurred  in  each  of  these 
instances,  it  is  assumed  that  the  resultant  shear  stress 
exceeded  the  shear  strength  given  by  the  Mohr-Coulomb 
failure  criterion  in  eq  5.  This  is  expressed  as 

rr  >  t0  +  on  tan0  (34) 


The  in  situ  stresses  on  the  in  situ  principal  planes 
are  represented  in  Figure  16b.  In  this  instance,  of 
course,  there  are  no  shear  components  and  the  prin¬ 
cipal  planes  lie  along  the  horizontal  axis.  Figure  16c 
represents  the  sum  of  the  dynamic  and  in  situ  hori¬ 
zontal  stresses  and  defines  a  new  set  of  principal 
stresses  resulting  from  a  combination  of  the  two.  The 
in  situ  principal  planes  lie  along  the  diameter  cd.  The 
normal  stress  on  the  major  principal  plane  is  now 
°rp  +  °pi>  wfl‘le  shear  component  remains  equal 
to  Tfp.  Likewise,  the  normal  component  on  the  minor 
plane  remains  r  ,  The  planes  upon  which  the  new 
combined  principal  stresses  act  are  oriented  at  the 
angle  w  with  respect  to  the  original  in  situ  principal 
planes.  Application  of  geometry  and  algebraic  man¬ 
ipulation  yields  an  expression  for  oj  of 


Also,  an  expression  for  the  major  and  minor  combined 
principal  stresses,  op  and  am  respectively,  can  be  de¬ 
rived: 


of  cos2o 


=  V"m,  +  0pl  A°r  + 

— ? - -It 

(gP! V1 *  +  2]'/l  (33) 


The  major  principal  stress  is  obtained  by  adding  the 
expression  in  brackets  and  the  minor  principal  stress 


Using  an  approximation  for  the  cohesion  obtained 
from  an  extrapolation  of  the  DIHEST  results,  the  max¬ 
imum  friction  angle  which  would  permit  motion  at 
each  slip  plane  location  on  MIGHTY  EPIC  was  calcu¬ 
lated. 

The  cohesion  approximation  is  derived  from  a  set 
of  assumptions  based  on  the  observation  (Blouin 
1980)  that  relative  displacements  were  never  observed 
beyond  a  range  of  three  crater  radii  during  the  DIHEST 
series.  Taking  this  obscrva*ion  as  a  generalization, 
the  cohesion  can  then  be  calculated  for  a  surface  burst 
geometry.  As  shown  in  the  schematic  diagram  of 
Figure  18,  the  peak  radial  stress  is  calculated  at  a 
range  of  three  crater  radii  from  a  surface  burst.  This 
is  assumed  to  be  the  maximum  range  at  which  slip 
will  occur.  As  in  the  case  of  the  contained  burst,  the 
critical  tangential  stress  is  taken  as  zero.  In  addition, 
it  is  assumed  that  near-surface  in  situ  stresses  are 
negligible.  Thus  the  peak  principal  stresses  of  zero 
and  a  ,  as  shown  on  the  Mohr  diagram  of  Figure  19, 
will  act  on  any  potential  slip  planes  at  that  range. 

Slip  will  occur,  as  shown,  if  a  slip  plane  is  oriented 
such  that  its  failure  envelope  is  just  tangent  to  the 
Mohr  circle.  This  condition  will  be  satisfied  only  by 
planes  oriented  at  an  angle  of  45°  ±0/2  to  the  radial 
stress  vector.  For  any  such  plane,  the  conditions  for 
tangcncy  to  the  circle  are  governed  by  a  combination 
of  0  and  tq.  From  geometry,  the  equation  for  co¬ 
hesion  as  aTunction  of  peak  radial  stress  and  friction 
angle  can  be  expressed  as 
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'  Mot.  Combined  S*r«»» 


Figure  1 7.  Schematic  plan  view  -combined  horizontal  principal  stresses. 


A-  L<p 

B-Shot  Point 
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EC3-Oi»ploc«d  Block 

Figure  18.  Schematic  section  view— hypothesized  limit  for  surface  burst 
relative  displacement. 


Figure  1 9.  Assumed  near-surface  critical  stress  stare. 
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The  value  of  radial  stress  at  a  range  of  three  crater 
radii  is,  of  course,  dependent  on  the  surface-burst 
crater  radius  which,  according  to  Crawford  et  al. 

(1974),  is  subject  to  a  great  deal  of  uncertainty,  par¬ 
ticularly  in  the  case  of  high  yield  nuclear  bursts.  This 
uncertainty  results  in  peculiar  variations  in  cohesion 
as  a  function  of  yield.  The  prediction  relationships 
recommended  by  Crawford  et  al.  (1974)  are  based  on 
the  high  yield  nuclear  tests  conducted  at  Eniwetok 
and  Bikini  atoll',  which  produced  abnormally  broad 
and  shallow  craters  as  compared  to  those  from  low 
yield  nuclear  and  HE  events  (Ristvet  et  al.  1978).  It 
is  not  known  whether  this  abnormality  was  geology- 
dependent  or  somehow  related  to  the  high  yields. 
However,  it  has  the  effect  of  making  the  cohesion 
according  to  eq  35  yield-dependent. 

Figure  20  is  a  plot  of  cohesion  as  a  function  of  yield 
for  high  energy  nuclear  and  HE  surface  bursts.  It  is 
based  on  the  prediction  techniques  for  dry,  soft  rock 
outlined  in  Crawford  et  al.  (1974)  and  in  eq  35.  The 
plot  includes  a  band  of  friction  angles  from  20  to  40°. 
This  encompasses  most  values  presented  in  the  litera¬ 
ture.  The  discontinuity  between  the  high  explosive 
and  nuclear  portions  of  the  plot  near  the  1-Kton  level 
is  due  to  several  factors.  First,  the  cratering  efficiency 
of  the  HE  is  about  20  times  that  of  the  nuclear,  and 
second,  the  equivalent  yield  coupling  factor  for  half- 
buried  HE  is  about  20%,  versus  only  4%  for  the  nuclear 
surface  burst.  Using  the  radial  stress  given  in  Figure  13, 
i.e.,  no  effect  of  the  free  surface  is  included  except 
the  redaction  in  coupling  efficiency  already  noted, 
the  band  of  cohesion  ranges  from  about  9  to  14  bars 
for  a  10-Mton  nuclear  surface  burst.  This  increases  to 
45  to  70  bars  for  a  1-Kton  nuclear  burst,  but  remains 
constant  with  yield  at  between  1 1  and  1 8  bars  for  the 
HE  bursts.  For  comparison,  the  calculated  cohesion 
from  the  ROCKTEST  ll-DIHEST  displacements 
listed  in  Table  1  ranged  from  35  to  53  bars,  which 
tends  to  agree  witn  the  low  yield  nuclear  portion  of  the 
plot.  It  should  be  remembered,  however,  that  all 
significant  DIHEST  related  displacements  occurred 
along  joints,  rather  than  faults,  which  would  tend  to 
have  a  higher  cohesion  because  the  small  interlocking 
asperities  have  not  been  sheared  off  by  previous  move¬ 
ment.  For  this  reason,  as  well  as  the  desire  to  keep 
the  assumed  MIGHTY  EPIC  cohesion  on  the  conser¬ 
vative  side,  a  cohesion  of  10  bars,  which  lies  on  the 
low  side  of  the  calculated  values,  was  used  in  the 
MIGHTY  EPIC  analysis  and  for  the  DIABLO  HAWK 
predictions. 

This  completes  the  engineering  approximations  and 
assumptions  necessary  to  evaluate  the  shear  strength 
parameters  al  each  displacement  location  on  MIGHTY 
EPIC.  In  summary,  the  Mohr-Coulomb  failure  crite¬ 
rion  is  used  to  obtain  the  upper  limit  of  the  friction 


angle  at  each  location  using  a  rearrangement  of  eq  34 
given  by 


^  <  tan'1 


(36) 


According  to  this  assumption,  the  in  situ  friction 
angle  must  have  been  less  than  or  equal  to  the  value 
calculated  at  each  location  for  slip  to  have  occurred. 

A  cohesion  of  10  bars  was  assumed,  based  on  the 
extrapolation  of  the  outlined  DIHEST  experience. 

The  values  of  the  peak  normal  and  resultant  shear 
stresses  acting  at  each  location  were  obtained  from 
eq  6  and  23-25.  The  peak  principal  stresses  used  in 
these  equations  were  calculated  by  superimposing  the 
critical  dynamic  principal  stresses  on  the  in  situ  stress 
approximation  of  Figure  8.  The  parameters  used  at 
each  MIGHTY  EPIC  displacement  location  arc  listed 
in  Table  3.  Because  they  arc  classified,  the  dynamic 
stresses  and  peak  combined  principal  stresses  are  not 
included  in  the  uble.  The  computed  uppe:  limits  of 
the  friction  angle  at  each  displacement  loc  don,  using 
the  assumed  value  for  cohesion  of  10  bars,  arc  listed. 
These  ranged  from  61.5°  at  the  3BC/3D  bedding 
plane  to  29.0°  at  the  intersection  of  fault  5  with  the 
Bypass  drift.  Again,  these  arc  the  highest  friction 
angle  values  which  would  have  allowed  slip  at  each 
location.  Had  the  values  been  less,  slip  would  still 
have  occurred.  Had  values  been  greater  than  those 
listed,  slip  would  not  have  occurred. 

The  direction  of  the  resultant  shear  vector  rf  is 
denoted  in  Table  3  by  the  angle  y.  This  angle  is 
measured  in  the  plane  of  the  fault,  viewing  the  fault 
from  the  working  point  side.  It  is  the  angle  between 
the  fault  strike  and  the  resultant  shear  vector  and  is 
given  by 

y  -  tan-1  (~^  (37) 

The  sign  convention  is  for  positive  values  of  y  to  have 
a  vertical  component  of  motion.  Left  lateral  motion 
means  that  the  fault  block  opposite  the  working 
point  block  should  move  to  the  left,  relative  to  the 
working  point  block.  Right  lateral  motion  infers  rel¬ 
ative  movement  to  the  right.  As  shown  in  Table  3, 
the  directions  of  motion  measured  on  MIGHTY  EPIC 
generally  agree  well  with  the  calculated  directions  of 
the  peak  shear  stress  vectors.  The  poorest  agreement 
is  the  38°  discrepancy  on  the  fault  which  intersects 
the  Interface  drift.  This  discrepancy  can  be  explained 
by  the  fact  that  the  working  point  block  must  have 
an  upward  component  of  motion  as  it  moves  outward 
along  the  underlying,  upward-sloping  quartzite  inter¬ 
face  shown  in  Figure  6.  Thus,  instead  of  reverse 
thrust  motion  on  the  interface  fault,  which  is  typical 
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Figure  20.  Cohesion  to  limit  relative  displacement  to  three  crater  radii. 
Table  3.  MIGHTY  EPIC  displacement  parameters. 
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of  all  oihei  MIGHTY  EPIC  and  DIHEST  data,  motion 
is  in  the  normal  mode. 

Displacement  analysis 

It  was  shown  in  the  DIHEST  analysis  section  that 
the  upper  bound  of  the  magnitude  of  displacements 
associated  with  the  DIHEST  shots  was  given  by  a 
conservation  of  block  energy  equation  (eq  4)  where 
the  block  has  an  initial  kinetic  energy  proportional 
to  the  square  of  the  peak  particle  velocity  at  the  range 
of  slip.  In  instances  where  geometric  or  air  overpressure 
constraints  did  not  inhibit  block  motion,  maximum 
differential  displacements  closely  matched  those  cal¬ 
culated  using  eq  4.  This  approach  was  modified  to 
account  for  deep  burial,  with  parameters  adjusted  to 
yield  good  agreement  between  the  MIGHTY  EPIC 
displacements  and  the  calculated  displacements.  This 
semi-empirical  technique  was  then  used  to  predict 
the  DIABLO  HAWK  block  motion  magnitudes. 

The  extrapolation  from  surface  relative  displace¬ 
ments  driven  by  near-surface  cratering  bursts  to  sub¬ 
surface  displacements  driven  by  contained  bursts  is 
complicated  by:  1 )  the  subsurface  displacement  blocks 
are  defined  only  at  their  intersection  with  post-test 
excavations,  and  2)  credible  ground  motion  and  stress 
time  histories  directly  adjacent  to  planes  of  slip  do 
not  exist.  The  first  complication  makes  it  impossible 


to  estimate  the  total  kinetic  energy  of  a  displaced 
block.  The  second  makes  it  difficult  to  assess  the 
normal  stress  distribution  on  a  fault  surface  during 
displacement,  which  is,  of  course,  the  key  to  deter¬ 
mining  the  frictional  energy  absorption. 

In  order  to  overcome  the  first  difficulty,  an  incre¬ 
ment  of  the  block  mass  is  used  which  acts  on  an  in¬ 
crement  of  the  fault  plane  as  shown  in  Figure  21 .  This 
mass  increment  is  assumed  to  extend  radially  from  the 
cavity  wall  formed  by  the  detonation  to  the  slip  plane. 
Its  massrrtj  is  given  by 


where  p  is  the  mass  density  of  the  rock,  the  radius 
of  the  increment  at  the  slip  plane,  r  the  range  from 
the  working  point  to  the  slip  plane,  and  rf  the  cavity 
radius  formed  by  the  explosion.  The  area  of  the  in¬ 
crement,  normal  to  the  radius,  from  the  working  point 
to  the  slip  plane,  given  by  iti  q,  is  taken  as  unity. 
Equation  3  can  now  be  rewritten  in  a  form  suitable  to 
the  underground  case  as 

f« 

Vtmy2m  •  pft  I  an  (x) dx  +  m,^6v  (39) 
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where  the  kinetic  block  energy,  V-jny*,  equals  half 
the  increment  mass  multiplied  by  the  square  of  the 
peak  radial  velocity  at  the  slip  plane.  The  first  group 
of  terms  to  the  right  of  the  equal  sign  represents  the 
energy  dissipated  by  friction.  This  is  the  integral  of 
the  normal  stress  an  (x)  acting  on  the  plane,  which 
is  assumed  to  vary  with  displacement  in  the  x  direc¬ 
tion  (taken  as  the  direction  of  the  principal  shear 
stress),  multiplied  by  the  coefficient  of  friction  p, 
and  a  term  f  which  equals  the  area  of  the  slip  plane 
intercepted  By  the  mass  increment.  The  area  f  is 
given  by 

*  sin?;  sin  0 

where  »j  is  the  angle  between  the  radial  velocity  vector 
and  the  slip  plane  strike  and  0,  the  slip  plane  dip.  Thus 
the  entire  integral  term  represents  a  varying  friction 
force  integrated  over  the  length  of  displacement  5. 

The  final  term  in  eq  39  is  simply  the  change  in  poten¬ 
tial  energy  due  to  vertical  translation  of  the  mass  in¬ 
crement.  By  defining  the  vertical  component  of  dis¬ 
placement  5y  as  a  function  of  the  total  displacement 
and  slip  plane  orientation,  and  by  approximating  the 
integral  terms,  the  total  displacement  6  can  be  obtained. 

By  approximating  the  integral  term  the  second 
complication  of  the  contained  burst  extrapolation, 
that  of  lack  of  appropriate  stress  time  histories,  is 
circumvented.  Three  different  approximations  of  the 
integral  term  were  used  to  make  post-test  "predic¬ 
tions’’  of  MIGHTY  EPIC  displacement  magnitudes. 

The  model  giving  the  best  overall  fit  to  the  displace¬ 


ment  data  was  then  used  to  make  the  DIABLO  HAWK 
predictions.  The  three  approximations  are  shown  in 
Figure  22.  In  all  cases  the  normal  stress  is  initially 
assumed  to  equal  the  peak  normal  stress  on  of  eq  25, 
computed  from  superimposing  the  peak  dynamic 
normal  stress  on  the  in  situ  normal  stress.  The  normal 
stress  at  the  conclusion  of  slip  is,  in  all  cases,  assumed 
equal  to  the  pretest  in  situ  normal  stress  onj  on  the 
slip  plane.  It  is  calculated  using  oq  25  with  opj  and 
omi  substituted  for  op  and  <Jm,  respectively,  and  with 
0  computed  as  the  angle  between  the  strike  of  the  slip 
plane  and  the  direction  of  opi.  The  differences  be¬ 
tween  the  three  models  are  in  the  rates  cf  normal 
stress  dissipation.  Model  1  assumes  that  the  dynamic 
component  is  dissipated  linearly  over  the  entire  dis¬ 
placement  6.  Model  2  assumes  that  the  dissipation 
occurs  linearly  in  only  half  the  total  displacement, 
with  the  remaining  displacement  occurring  under  the 
in  situ  normal  stress.  Model  3  assumes  a  linear  dissi¬ 
pation  which  is  independent  of  the  total  displacement. 
Trial  and  error  showed  that  a  displacement  of  0.3  m 
gave  the  most  reasonable  fit  to  the  MIGHTY  EPIC 
data.  Again,  the  remaining  slip  was  assumed  to  occur 
ur  der  the  in  situ  normal  stress.  The  integral  in  eq  39 
was  thus  defined  as  the  area  within  each  of  the  three 
normal  stress  distributions  shown  in  Figure  22.  It 
should  be  noted  that  none  of  the  three  models  is 
meant  to  rigorously  represent  the  actual  normal  stress 
attenuation  on  the  plane  of  weakness.  They  arc  simple 
representations  of  possible  attenuation  mechanisms 
and  provide  a  simple  means  of  evaluating  the  integral 
in  eq  39. 

The  vertical  component  of  relative  displacement 
8V  is  given  by 


Working 

Point 


Figure  21.  Depiction  of  mass  increment. 


'formal  stress  dissipation  models. 
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Win0  (41) 

where  5rf  is  the  dip  component  of  displacement. 
Assuming  displacement  occurs  in  the  direction  of  the 
peak  resultant  shear  stress,  Sd  can  then  be  defined  as 

5d=5sin7  (42) 

where  y  is  the  angle  between  the  resultant  shear  stress 
and  the  slip  plane  strike  shown  in  Figure  1 2.  By  inspec¬ 
tion  •  • 

sin  y=  Il_  (43) 

Combining  eq  41-43  gives  an  expression  for  the  verti¬ 
cal  component  of  displacement  as  a  function  of  the 
dip  shear  stress,  resultant  shear  stress  and  slip  plane  dip: 


responds  to  a  friction  angle  of  29°,  which  is  the  lowest 
limit  computed  in  the  shear  stress  analysis  of  the 
MIGHTY  EPIC  relative  displacements.  The  computed 
area  coefficients,  in  situ  normal  stress  and  shear  stress 
ratios  are  listed  in  Table  4,  but,  because  of  classifica¬ 
tion  restrictions,  the  values  of  mass  increments,  peak 
normal  stresses  and  peak  velocities  are  omitted.  The 
displacements  computed  using  each  of  the  three  models 
are  listed  along  with  the  actual  displacements.  Overall, 
model  2  (eq  46),  which  assumes  that  the  dynamic 
component  of  normal  stress  dissipates  over  half  the 
total  displacement,  gives  the  best  agreement  with  the 
measurements.  All  values  of  computed  displacement 
using  this  method  fall  within  about  40%  of  the  actual 
values. 


DIABLO  HAWK  PREDICTIONS 


Substituting  the  areas  of  each  of  the  normal  stress 
dissipation  models  for  the  integral  in  eq  39  and  sub¬ 
stituting  eq  44  for  6y  gives  solutions  of  total  displace¬ 
ment  for  each  of  the  dissipation  models.  The  solution 
for  model  1  in  which  the  dynamic  normal  stress  is 
linearly  dissipated  over  the  entire  displacement  is 
given  by 

5  = - j -  (45) 

^K+0ni)+2"V7?pSin/J 


The  solution  for  model  2  in  which  dynamic  normal 
stress  is  dissipated  over  half  the  total  displacement . 
is  given  by 

m v* 

5= _ LI _  (46) 

Vtufu  (on  +  3  an. )  +  2/^—1  sind? 

The  solution  for  model  3  in  which  dynamic  normal 
stress  is  dissipated  over  0.3  m  displacement  is  given 
by  ' 


„  Vvntf  '  °-15  Vfg(on  '  °n| ) 

J 

Uf%  +  ™  sin/J 


(47) 


Each  of  these  three  displacement  equations  was 
applied  to  the  five  documented  relative  displacements 
from  MIGHTY  EPIC  listed  in  Table  4,  A  coefficient 
of  friction  (i  of  0.55  was  used  throughout.  This  cor¬ 


The  DIABLO  HAWK  relative  displacement  predic¬ 
tions  are  based  on  the  MIGHTY  EPIC  analysis  out¬ 
lined  in  the  previous  section.  Like  that  analysis,  it  is 
divided  into  two  parts.  The  first  part  includes  a  shear 
stress  analysis  which  computes  the  anticipated  magni¬ 
tudes  and  directions  of  the  peak  shear  stresses  at  each 
DIABLO  HAWK  fault  documentation  location  and 
compares  these  to  the  predicted  shear  strengths.  At 
locations  where  the  predicted  shear  stress  exceeds  the 
shear  strength,  relative  displacement  is  predicted  in 
the  direction  of  the  peak  shear  stress.  The  second 
portion  of  the  analysis  is  then  applied  only  at  these 
locations  to  predict  the  magnitude  of  the  relative 
displacements. 

Figure  23  is  a  plan  view  of  all  passive  fault  docu¬ 
mentation  locations  on  DIABLO  HAWK.  The  faults 
are  identified  by  the  drift  in  which  they  are  located, 
according  to  the  system  used  by  Kipp  and  Kennedy 
(1978).  The  strike,  dip,  radial  distance  and  direction 
from  the  working  point  of  the  faults  at  each  docu¬ 
mentation  station  are  listed  in  Table  5.  In  instances 
where  the  faults  can  be  identified  as  part  of  a  major 
throughgoing  fault,  they  are  also  listed  by  the  number 
of  the  major  fault  adjacent  to  the  drift  identification 
number. 

The  passive  fault  documentation,  prepared  by  the 
U.S.  Geological  Survey,  consists  of  four  steel  pins 
glued  in  the  form  of  a  parallelogram  into  the  drift 
wall  as  shown  schematically  in  Figure  24.  Two  pins 
are  located  on  each  side  of  the  fault,  roughly  paral¬ 
lel  to  it.  This  arrangement  is  duplicated  on  the 
opposite  wall  of  the  drift,  making  a  total  of  eight 
pins  at  each  fault  intersection.  The  28  different  pin 
to  pin  dimensions  are  recorded  before  and  after  the 
test.  Triangulation  is  used  to  determine  the  strike 
and  dip  components  of  any  relative  displacement. 
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Figure  24.  Schematic  side  viev.-lypicat 
pasuve  tauit  documentation  station 


The  determination  of  the  peak  shear  stresses  jt 
each  fault  documentation  location  follows  the  tech¬ 
nique  outlined  in  the  presious  section.  The  predicted 
peak  values  of  strike  shear,  dip  shear  and  normal 
stress  were  determined  using  eq  23,  24  and  25  rcspec- 
r.elv  .  The  peats  principal  stresses  used  in  these  equa¬ 
tions  were  computed  hy  superimposing  the  critical 
dynamic  stresses  on  the  in  situ  stress  approximation 
tor  Rainier  Mesa.  As  in  the  MIGHTY  T P I C  analysis, 
the  critical  dynamic  stress  slate  was  assumed  to  be  the 
peak  radial  stress  given  m  Figure  13  acting  simulta¬ 
neously  with  ?cro  tangential  stress.  The  superposition 
was  done  using  Mohr  circles,  In  addition  to  the  peak 
principal  stresses,  the  Mohr  circle  analysis  also  y  ielded 
values  fur  the  angle  0  required  (or  the  solution  of  eq 
23-25  at  each  location.  The  peak  strike  ana  dip  shear 
stresses  were  substituted  into  cq  6  to  give  the  resultant 
total  shear  stress  acting  in  the  direction  7  as  listed  in 
Table  5.  The  peak  normal  stresses  were  substituted 
into  eq  5  to  obtain  the  critical  shear  stresses  (antici¬ 
pated  shear  strengths)  at  each,  location.  The  parameters 
obtained  from  the  DIHEST  and  MIGHTY  C P  1C 
analyses  (i.c.,  a  cohesion  r0  of  1 0  bars  and  a  friction 
angle  of  29°!  were  used  to  calculate  7  .  The  rationale 
in  employing  these  values,  discussed  in  detail  in  the 
previous  section,  was  that  1 0  bars  appeared  to  approx¬ 
imate  the  lower  bound  of  an  extrapolation  of  the 
DIHTST  lesulls  to  nui  iear  and  high  explosive  surface 
bursts,  and  the  friction  angle  of  29  Is  the  computed 
maximum  value  which  would  have  permitted  all  the 
displacements  observed  on  MIGHTY  EPIC  to  occur. 

The  ratios  of  the  peak  shear  stress  to  the  critical 
shear  strew,  7  r  ,  are  lifted  in  Table  5.  f  or  values 
greater  than  1.0,  the  total  shear  stress  exceeds  the 
computed  shear  strength  and  relative  displacement  is 
predicted.  The  stress  ratio  may  be  indicative  of  the 
probability  of  displacement  occurrence.  Higher  ratios 


would  indicate  a  higher  likelihood  of  displacement 
than  ratios  closer  to  1 .0.  Of  the  laulls  lhal  arc  not 
predicted  to  move,  those  with  a  ratio  close  to  1 .0 
would  have  a  higher  likelihood  of  displacing  than  those 
with  lower  stress  ratios.  Again,  because  of  classifica¬ 
tion  restrictions,  the  actual  values  of  shear  stresses 
and  normal  stresses  are  not  included  in  the  table. 

Displacement  predictions  were  made  for  the  14 
fault  documentation  stations  with  stress  ratios  greater 
than  1 .0  using  cq  46,  which  gave  the  best  fit  to  the 
MIGHTY  EPIC  data  of  the  three  displacement  models 
examined.  The  computed  area  coefficients  f '  .  in  situ 
normal  stresses  0  and  dip  shear  stress  to  total  stress 

n  1 

ratios  r./Tf  are  listed  in  Tabic  S.  The  remaining  in¬ 
put  parameters  are  not  listed  because  of  classification 
restrictions  The  14  total  displacement  predictions 
are  listed  along  with  a  breakdown  of  their  strike  and 
dip  components.  The  strike  and  dip  components  are 
compuied  with  the  assumption  that  motion  will  be  in 
the  direction  of  the  resultant  shejr  vector  7  ■  Negative 
dtp  slip  values  indicate  that  the  block  on  the  working 
point  side  of  the  fault  moves  downward  relative  to 
the  opposite  block.  RL  and  LL  indicate  right  lateral 
or  left  lateral  strike  slip  displacement.  The  predicted 
magnitude  of  total  displacement  ranges  Irom  0.1  5  in 
at  BP-4  to  1 .83  m  at  BP-9.  The  shear  stress  ratios 
for  several  ol  the  faults  tor  which  displacement  is  pre¬ 
dicted,  such  as  BP — 4  and  RE-1 ,  are  close  enough  to 
5.0  that  it  would  not  be  surprising  if  no  relative  mulic  n 
uccurrcd.  On  the  other  hand,  there  are  several  faults, 
such  as  LOS-  3  and  AB-6,  for  which  no  displacement 
is  indicated.  But  they  have  stress  ratios  close  to  1 .0 
and  might  easily  exhibit  relative  displacement. 

In  summary,  the  DIABLO  HAWK  predictions  are 
based  on  semi-empirical  equations  derived  from  the 
combined  MIGHTY  EPIC  and  DIHEST  experience. 

The  estimated  peak  shear  stress  acting  on  each  fault 
is  compared  to  the  estimated  shear  strength  to  deter¬ 
mine  whether  slip  will  occur  Both  estimates  arc- 
based  on  approximations  of  the  combined  in  situ 
stress  and  dynamic  loadings  If  slip  is  predicted,  the 
magnitude  of  displacement  is  computed  from  an 
equation  based  on  the  dissipation  ol  the  kinetic  energy 
for  an  increment  of  mass  extending  radially  from  the 
cavity  wall  to  the  fault.  The  equation  considers  energy 
dissipated  by  friction  along  the  fault  duiing  displace¬ 
ment  and  the  loss  or  gain  in  energy  associated  with 
changes  in  elevation  of  the  block.  Both  the  coeffi¬ 
cient  of  sliding  friction  and  the  norma!  stress  distri¬ 
bution  on  the  fault  during  the  course  ol  the  relative 
motion  are  based  on  analyses  of  the  MIGH  I  Y  EPIC, 
data. 

Obviously  ,  the  principal  inadequacy  of  these  pre¬ 
diction  techniques  is  lhal  they  do  not  treat  the  faults 
and  blocks  making  up  the  in  situ  mass  as  continuous. 
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Table  5.  DIABLOW  HAWK  displacement  predictions. 


Predicted 
displacement  (m) 


Documentation 

station 

LOS- 1 
LOS-2 
LOS-3 
LOS-4 
BP-1 
BP-2 
BP-3 
BP -4 
BPS 
BP-6 
BP-7 
BP -8 
BP-8A 
BP -9 
BP-10 
RE-1 
AB-1 
AB-2 
AB-3 
AB-4 
AB-5 
AB-6 
CB-1 
CB-2 
CB-3 
CB-4 
CBS 


Malar 
fault  no. 


NI8W 

NSOE 

N28E 

N43W 

N25W 

N9W 

N2E 

N34W 

N29W 

N11E 

NS  IE 

N19E 

N45W 

N34W 

N20E 

N39E 

N5E 

N29W 

N24W 

N40W 

N21W 

N10W 

N24E 

N8E 

N2W 

N34E 

N11W 

N5E 


63SW 

82SE 

6SSE 

8SNE 

60NE 

7SNE 

65  SE 

85SW 

67SW 

75NW 

80SE 

70NW 

80NE 

80SW 

60SE 

79SE 

60SE 

65SW 

70SW 

7SSW 

7SNE 

65NE 

68  NW 

80NW 

72NE 

70SE 

67SW 

70NW 


Direction 
from  W.P. 

rr/Tc 

> 

(degl 

'i 

°n, 

(bars) 

Td/rr 

*< 

strike • 

kd 

dip 

6 

Total 

S79E 

0.94 

S79E 

1.21 

-  8.6 

1.35 

40.3 

-0.15 

0.5 1  RL 

-0.08 

0.52 

S79E 

0.90 

S79E 

2.06 

-  3.4 

1.67 

80.0 

-0.06 

1.07  LL 

-0.06 

1.07 

S72E 

0.54 

S72E 

0.41 

S71E 

0.46 

S70E 

1.04 

4.0 

1.56 

78.6 

0.07 

0.15  LL 

0.01 

0.15 

S70E 

1.16 

19.7 

1.66 

77.4 

0.33 

0.14  LL 

0.05 

0.15 

S69E 

0.40 

S66E 

0.86 

S6SE 

0.57 

S66E 

2.47 

-  3.4 

2.71 

80.0 

-0.06 

0.54  LL 

-0.03 

0.54 

S56E 

3.59 

4.0 

2.77 

78.6 

0.07 

1.82  LL 

0.13 

1.82 

S76W 

1.68 

35.5 

1.40 

SS.i 

0.58 

1.29  RL 

0.92 

1.58 

N84W 

1.19 

17.5 

’..23 

42.2 

0.30 

1.20  RL 

0.35 

1.15 

S60W 

1.62 

35.5 

1.40 

62.4 

0.58 

0.96  RL 

0.68 

1.18 

1.21  50.0 

8.97  55.3 

2.10  63.2 


0.59  1.21  LL  0.88  1.50 

0.0  1.36  LL  0.0  1.36 

0.14  0.55  RL  0.08  0.56 

-0.33  0.33  RL  -0.12  0.35 

-0.15  0.34  RL  -O.OStU* 


*  RL  is  right  lateral  strike  slip  displacement;  LL  is  left  lateral  displacement. 


For  instance,  in  Figure  23  the  faults  labeled  BP-10, 

AB-1  and  CB-4  arc  all  segments  of  the  major  fault 
5.  At  documentation  stations  BP-10  and  AB-i  sub¬ 
stantial  slip  is  predicted.  Yet  at  station  CB-4,  because 
the  fault  s'rikc  has  changed  so  that  it  trends  directly 
toward  the  working  point,  motion  is  not  predicted. 

In  fact,  the  shear  stress  ratio  is  very  low,  which  would 
normally  be  taken  to  mean  that  the  probability  of 
displacement  is  also  very  low.  However,  if  substan¬ 
tial  displacement  occurs  as  predicted  at  both  BP-10 
and  AB-1  it  would  seem  likely  that  this  would  carry 
through  to  CB-4  unless  some  other  intervening  fault 
provided  a  more  convenient  release  for  this  motion. 
Likewise,  the  magnitude  of  displacement  at  one 
point  within  a  displaced  fault  block  cannot  be  inde¬ 
pendent  of  the  magnitude  at  another  point  within 
that  same  block,  if  the  block  remains  intact  between 
the  two  points.  In  other  words,  if  a  meter  of  slip 
were  to  occur  at  BP-10  and  AB-1  due  to  displace¬ 
ment  of  the  block  on  the  working  point  side  of 
fault  5  a  similar  amount  of  displacement  would  be 
exoected  at  CB-4,  assuming  that  the  slip  block  is  an 
intact  block  and  that  residual  compression  of  the  block 
is  negligible. 


Unfortunately,  il  appears  impossible  to  adequat.  -  ■ 
describe  the  complete  geometry  of  the  planes  of 
weakness  and  blocks  surrounding  a  detonation,  let 
alone  the  intricate  details  of  their  interactions  during 
relative  displacement.  In  many  instances  a  potential 
slip  plane  would  only  be  encountered  at  one  point  by 
an  excavation.  The  directions  of  the  extension  of 
this  plane  and  the  existence  of  other  interacting 
planes  could  only  be  determined  through  costly  and 
elaborate  exploration  programs.  Even  with  the  most 
complete  exploration  programs  there  arc  often  struc¬ 
tural  weaknesses  of  critical  importance  that  go  unde¬ 
tected.  In  other  words,  it  appears  unlikely  that  a 
reliable  quantitative  model  which  accurately  de¬ 
scribes  detailed  block  motion  behavior  can  be  de¬ 
veloped  in  the  near  future  because  the  analysis  tools 
are  inadequate  for  the  three-dimensional  job  and 
the  geologic  exploration  necessary  to  completely 
describe  a  site  would  be  prohibitively  expensive. 
Perhaps  a  semi-empirical  approach,  coupled  with  a 
probabilistic  analysis  based  on  our  expanding  data 
base,  will  prove  adequate. 
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